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ABSTRACT.— The Asian newt genus Tylototriton, is the only urodelan genus currently recorded in Thailand. For 
a long time in the past, T. verrucosus was the sole known urodelan species found in Thailand, but was recently 
shown to be comprised of three cryptic species, T. uyenoi, T. anguliceps, and T. panhai. However, in our 
molecular study, the newts from Doi Chang, Chiang Rai Province were different from these three recognized 
Tylototriton species, and nested within the clade of T. verrucosus from China. They also had morphological traits 
similar to T. verrucosus. These results lead to the confirmation of the existence of 7. verrucosus in Thailand. The 
habitats of T. verrucosus need to be protected by future conservation plans. 


KEY WORDS: conservation, cryptic species, newt, Thailand, Tylototriton verrucosus 


INTRODUCTION 


Thailand is a rather small country that 
has undergone extensive land development, 
but is still relatively rich in biodiversity 
(Baimai, 2010). Several new species of 
anurans have been discovered in 
unprotected areas (e.g. Suwannapoom et al., 
2016, 2018). However, the majority of 
hotspot areas are inside government 
protected enclosures (Tantipisanuh and 
Gale, 2018), including national parks (NP) 
or wildlife sanctuaries (WS). Some of these 
protected areas house the newt genus 
Tylototriton, which have a small range 
distribution about which there is scarce 
information. 

Tylototriton, is the only genus of 
salamandrid, Urodela, currently recorded in 
Thailand. In the past, T. verrucosus was the 
sole known representative species found in 
Thailand (Taylor, 1962; Wongratana, 1984; 


Nabhitabhata et al., 2000; Chan-ard, 2003). 
Their habitats are ponds that are covered 
with aquatic plants along with sand, stone 
debris, and organic materials scattered on 
the floor. They also reside terrestrially in 
moist areas, such as under rock piles and 
piles of dead leaves. They are distributed 
along high mountain ranges in the North and 
North-eastern regions of Thailand at an 
altitude of more than 1000 meters above 
medium sea level (m amsl) (Taylor, 1962; 
Wongratana, 1984; Chan-ard et al., 1999, 
Chan-ard, 2003; Pomchote et al., 2008; 
Hernandez et al., 2019). From 1962 to 2007, 
all Tylototriton found in Thailand were 
classified as T. verrucosus, although they 
were clearly distinct in terms of 
geographical separation and color patterns. 
Pomchote et al. (2008) reported two types of 
T. verrucosus: Type I had a yellow to 
orange body coloration and Type II had a 
dull body coloration. They were found along 
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mountain ranges in the Northern and North- 
eastern parts of Thailand, respectively. 

The above research led to the 
subsequent studies that reported in fact T. 
verrucosus in Thailand was actually 
comprised of three species of newts: T. 
uyenoi, T. anguliceps, and T. panhai. These 
three species were distributed along the 
Ang Khang, Chang Kien, Inthanon, Pui, 
and Suthep mountains in Chiang Mai 
Province (T. uyenoi) (Nishikawa et al., 
2013a), the Lahnga mountain in Chiang 
Rai Province (T. anguliceps) (Le et al., 
2015), and the Phu Hin Rong Kla NP in 
Phitsanulok Province, and Phu Luang WS 
and Phu Suan Sai NP in Loei Province (T. 
panhai) (Nishikawa et al., 2013a). Despite 
confirmation of the existence of three newt 
species in Thailand, the question remained 
as to the existence of the true T. verrucosus. 
Should its existence be confirmed, where is 
it distributed? This question arose since, for 
example, a study conducted in 2013 
assigned T. uyenoi to T. verrucosus (in 
Taylor, 1962) and T. panhai had been 
reported as T. verrucosus (in Wongratana, 
1984; Chan-ard, 2003; Khonsue et al., 
2010). This led to the emphasis of the 
question whether or not T. verrucosus does 
exist in Thailand, or simply represents 
misidentifications since 1962 of T. uyenoi, 
T. panhai, and T. anguliceps with no T. 
verrucosus in Thailand. 

From our field survey in the Northern 
parts of Thailand, we found newts from Doi 
Chang, Chiang Rai Province, which is a new 
locality of Tylototriton in Thailand. 
However, this location is not in the 
protected area where the previous three 
Tylototriton species were observed. The 
newts found in Doi Chang had a rather 
similar color pattern to 7. uyenoi and T. 
anguliceps, but were rather different in other 
external characters. Furthermore, Doi Chang 


(Chiang Rai Province) is separated from the 
locations where T. uyenoi and T. anguliceps 
were previously found. From the above lines 
of information, the possibility of a different 
Tylototriton species existing in Doi Chang, 
Chiang Rai Province, is suggested. 
Therefore, molecular and morphological 
analyses were performed to examine the 
taxonomic status of Tylototriton species 
from Doi Chang. The knowledge obtained 
in this study will be used in the 
implementation of future conservation plans 
of Thai Tylototriton. 


MATERIALS AND METHODS 


sampling 

Field surveys using the visual encounter 
survey method (Heyer et al., 1994) were 
conducted from the 23" to the 24^ and the 
31* May 2018 in Namtok Mae Surin NP, 
Mae Hong Son Province and Khun Chae 
NP, Chiang Rai Province, respectively, and 
from the 2™ to the 3% June 2018 in Doi 
Chang, Chiang Rai Province (Fig. 1). Ponds 
selected as survey sites were those that had 
bottoms covered with small pieces of debris 
on the mountains with an elevation of at 
least 1000 m amsl (Pomchote et al., 2008). 
Upon an encounter, adult newts were 
captured by either an aquatic dip net or by 
hand and kept in plastic boxes. The 
biological and physical data of their habitats 
were also recorded. 

Male newts had obvious sexual cloacal 
characters. During the breeding season, the 
male cloaca is swollen and the cloacal 
opening is oval shaped (Pomchote et al., 
2008; Phimmachak et al., 2015b). From 
these criteria, all the newts found in this 
study were sexed as breeding males. Twenty 
newts (three from Khun Chae NP, seven 
from Doi Chang, and 10 from Namtok Mae 
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FIGURE 1. Distribution map of the three Tylototriton species in this study: Tylototriton anguliceps (circle), 
T. uyenoi (square), and T. verrucosus (triangle). The map is modified by Praphatsorn Punsompong. 


Surin NP) were used for the morphological The samples were anesthetized by 
study, while 15 newts were used for the immersion in a solution of tricaine methane 
molecular analyses (Table 1). sulfonate (MS-222; 5 g/L) for about five 


minutes (modified from Mitchell, 2009). All 
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TABLE 1. Specimens for molecular analyses. CAS = California Academy of Sciences; CIB = Chengdu Institute of 
Biology; HNUE = Hanoi National University of Education; IEBR = Institute of Ecology and Biological Resources, 
Hanoi; KIZ = Kunming Institute of Zoology; KUHE = Graduate School of Human and Environmental Studies, Kyoto 
University; LSUHC = La Sierra University Herpetological Collection; MVZ = Museum of Vertebrate Zoology, 
University of California, Berkeley; NMNS = National Museum of Natural Science, Taiwan; TBU = Tay Bac 
University, Son La; VNMN = Vietnam National Museum of Nature; ZDUM = Zoology Department of University of 
Mandalay; ZMMU = Zoological Museum of Moscow State University. *topotype 


Sample Species Voucher no. Locality GenBank Source 
no. accession no. 
1 T. verrucosus * KIZ Husa, Yunnan, China AB922818 Nishikawa et al. 
201306055 (2014) 
2 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505605 This study 
Thailand 
3 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505606 This study 
Thailand 
4 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505607 This study 
Thailand 
5 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505608 This study 
Thailand 
6 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505609 This study 
Thailand 
7 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505610 This study 
Thailand 
8 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505615 This study 
Thailand 
9 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505616 This study 
Thailand 
10 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505614 This study 
Thailand 
11 T. verrucosus No voucher Doi Chang, Chiang Rai, LC505617 This study 
Thailand 
12 T. anguliceps* VNMN Muong Nhe, Dien Bien, LC017832 Le et al. (2015) 
A.2014.3 Vietnam 
13 T. nguliceps* HNUE Muong Nhe, Dien Bien, LC017836 Le et al. (2015) 
A.I.1.111 Vietnam 
14 T. anguliceps TBU Thuan Chau, Son La, LC017833 Le et al. (2015) 
PAE.671 Vietnam 
15 T. anguliceps IEBR Thuan Chau, Son La, LC017834 Le et al. (2015) 
A.2014.10 Vietnam 
16 T. anguliceps No voucher Doi Lahnga, Chiang Rai, AB830728 Nishikawa et al. 
Thailand (2013a) 
17 T. anguliceps No voucher Khun Chae NP, Chiang Rai, LC505611 This study 
Thailand 
18 T. anguliceps No voucher Khun Chae NP, Chiang Rai, LC505613 This study 
Thailand 
19 T. uyenoi * KUHE 19146 Doi Suthep, Chiang Mai, AB830732 Nishikawa et al. 
Thailand (2013a) 
20 T. uyenoi * KUHE 19147 Doi Suthep, Chiang Mai, AB830733 Nishikawa et al. 
Thailand (2013a) 
21 T. uyenoi * No voucher Doi Suthep, Chiang Mai, AB830734 Nishikawa et al. 
Thailand (2013a) 
22 T. uyenoi KUHE 19037 Doi Inthanon, Chiang Mai, AB830730 Nishikawa et al. 


Thailand 


(20132) 
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procedures were then conducted once the 
sample was fully anesthetized. After 
completing the morphological measurements 
and tissue sampling for molecular analysis 
(see below), all anesthetized newts were 
rinsed in water to help them regain 
consciousness. They were assessed by 
observation of their movements, kept and 
monitored in a captive box for 24 h to allow 
full recovery before being released on land 
near their respective capture site. 

Molecular Study 

Tissue samples were obtained by 
removing the tail tips, approximately 8% of 
the full tail length (Arntzen et al., 1999; 
Polich et al., 2013), using surgical scissors 
sterilized with 70% (v/v) ethanol (Fellers 
and Freel, 1995). Then an antiseptic liquid 
was used to clean the wound. The tissue 
samples were preserved in absolute ethanol 
for DNA extraction and further analysis. 

Total DNA was extracted using a DNeasy 
Blood & Tissue kit (QIAGEN, Hilden, 
Germany) and used to PCR amplify and 
sequence the mitochondrial NADH 
dehydrogenase subunit 2 region (ND2) for 
molecular analyses as reported previously 
(Nishikawa et al., 2013a) Sequence data 
obtained were deposited in GenBank 
(accession Nos. LC505605—619). Combining 
them with sequences of the other Tylototriton 
species (Table 1), molecular phylogenetic 
trees were reconstructed by using the 
maximum likelihood (ML) and Bayesian 
inference (BI) methods. We used Kakusan 4 
(Tanabe, 2011), RAxML v8.2.8 (Stamatakis, 
2014), and MrBayes v3.2.6 (Ronquist et al., 
2012) for model selection, ML, and BI, 
respectively. 

The robustness of the ML tree was tested 
using bootstrap analysis (Felsenstein, 1985) 
with 1000 replicates, and we regarded tree 
topologies with bootstrap values (BS) of 
70% or greater as well supported. For BI, 





two independent runs of four Markov chains 
were conducted for 10 million generations. 
The first three million generations were 
discarded as burn-in, and we considered 
topologies with posterior probabilities (bpp) 
of 0.95 or greater as significant support. 
Pairwise comparisons of uncorrected 
sequence divergences (p-distance by 400 
bps) were calculated using MEGA X 
(Kumar et al., 2018). 
Morphological Study 

The morphology of the newts was 
compared among the three locations (total N 
= 20). The following 27 measurements, with 
character definitions are described in 
Nishikawa et al. (2011), were taken for 
comparisons: SVL (snout-vent length); HL 
(head length); HW (head width); MXHW 
(maximum head width); SL (snout length); 
LJL (lower jaw length); ENL (eyelid-nostril 
length); IND (internarial distance); IOD 
(interorbital distance); UEW (upper eyelid 
width); UEL (upper eyelid length); OL 
(orbit length); AGD (axilla-groin distance); 
TRL (trunk length); TAL (tail length); VL 
(vent length); BTAW (basal tail width); 
MTAW (medial tail width); BTAH (basal 
tail height); MXTAH (maximum tail 
height); MTAH (medial tail height); FLL 
(forelimb length); HLL (hindlimb length); 
2FL (second finger length); 3FL (third 
finger length); 3TL (third toe length); and 
STL (fifth toe length). All measurements 
were taken using a digital sliding caliper to 
the nearest 0.01 mm. Each measurement 
was conducted three times and the average 
was used for further analyses. While the 
newts were anesthetized, their body weights 
(BW) were recorded using a digital 
weighing scale. The sex was checked from 
observing the cloacal characters (noted in 
Results). 

The SVL, BW, and 26 ratio values to 
SVL (R, %) were compared among the 
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newts from the three locations. Differences 
in morphological characters between 7. 
anguliceps and T. verrucosus, and between 
T. uyenoi and T. verrucosus were analyzed 


by the Mann-Whitney U test. The 
relationships of all characters were 
61 
99 

97 

69 
90 
64 


examined using principal | component 
analysis (PCA). All statistical analyses were 
performed using SPSS version 22 for 
Windows. Statistical significance was 
accepted at the P « 0.05 level. 
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FIGURE 2. Maximum likelihood phylogenetic tree based on the partial ND2 mtDNA gene for Tylototriton 
and related species (see Table 1). Numbers at nodes show bootstrap supports for ML inference (bs). 
Numbers at branches tips are sample numbers shown in Table 1. Scale bar — 0.02 substitutions/site. 
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RESULTS 


Molecular Study 

We analyzed ~1035 bp sequences of the 
partial ND2 region for 49 specimens, 
including an outgroup from GenBank data 
(Table 1). Of the 1035 nucleotide sites 
sequenced in this study, 369 were variable 
and 262 were parsimony informative. The 
mean likelihood score of the BI analyses for 
all trees sampled at stationarity was - 
6078.038. The likelihood value of the ML 
tree was -6074. 

Phylogenetic analyses employing two 
different optimality criteria (BI and ML) 
yielded nearly identical topologies. We, 
therefore, present only the ML tree in Figure 2. 
However, whilst monophyly of Tylototriton 





was supported in the BI tree (bpp = 0.96), it 
was not supported in the ML tree (BS = 
69%). Within Tylototriton, we recognized 
many clades but their relationships did not 
tend to be well supported (genetic distances 
available upon request from the senior 
author). Thai samples from Doi Chang 
(Samples 2—11) were nested within Clade A 
(see Fig. 2) of T. verrucosus and T. shanjing 
(Samples 1 and 31) with an average genetic 
p-distance of 1% (range: 0—2%), those from 
Khun Chae NP (Samples 17 and 18) were 
nested within Clade B (see Fig. 2) of 7. 
anguliceps (Samples 12—16) with an average 
p-distance of 1% (range: 0-2%), and those 
from Namtok Mae Surin NP (Samples 24— 
26) were nested within Clade C (see Fig. 2) 
of T. uyenoi (Samples 19—23) with an 











TABLE 2. Morphological comparison of adult males of the three Tylototriton species from this study (mean + 
SD of SVL (in mm), mean + SD of BW (g), and median of ratios of characters (R: %SVL), with range in 


parentheses). Character abbreviations refer to text. 


Species T. anguliceps T. verrucosus T. uyenoi 
No. and character 3 7 10 
SVL 64.8 + 3.5 66.6 x 2.7* 71.85 2.4 
(60.9-67.4) (62.6-70.3) (68.9-75.8) 
BW 14.0 x 2.3 10.7 € 0.6* 14.82 1.8 
(12.2-16.6) (10.2-12.0) (11.2-17.0) 
RHL 26.1 26.4* 25.1 
(25.7-26.4) (24.9-27.5) (24.2-26.3) 
RHW 19.6 17.6* 18.8 
(19.5-20.3) (17.5-18.4) (17.5-19.3) 
RMXHW 24.4 24.17 25.8 
(24.4) (23.2-25.1) (24.5-26.4) 
RSL 74 8.17 8.8 
(7.1-8.2) (7.6-9.0) (8.1-9.4) 
RLJL 20.5 20.7* 22.0 
(20.1-21.4) (19.7-21.5) (20.7-22.5) 
RENL 6.2 6.2 6.8 
(6.1-6.3) (5.4-7.9) (6.0-7.5) 
RIND 7.8 6.2 6.8 
(6.8-7.9) (5.3-6.7) (5.6-7.5) 
RIOD 12.3 13.3 13.0 
(12.1-12.4) (12.1-13.5) (12.6-14.4) 
RUEW 32 2.6* 3.1 
(2.7-3.4) (2.2-3.1) (2.2-3.8) 
RUEL 7.3 6.2 6.4 
(6.8-7.3) (5.6-7.4) (5.8-7.1) 
ROL 4.2 4.4 4.2 
(4.2-4.4) (3.6-5.7) (3.5-4.8) 
RAGD 52.3 55.]* 49.9 
(51.1-52.4) (53.9-56.4) (45.7-52.3) 


"P « 0.05 compared to T. anguliceps (Mann-Whitney U test) 


*P « 0.05 compared to 7. uyenoi (Mann-Whitney U test) 





Species T. anguliceps T. verrucosus T. uyenoi 

No. and character 3 7 10 

: RTRL 712 73.67 75.0 
: (70.8-71.5) (72.1-74.5) (71.8-98.0) 

i RTAL 873 94.3" 98.0 
(87.2-87.8) (93.5-110.4) (88.8-110.4) 

i RVL 114 10.1 124 
(9.6-11.4) (9.7-11.9) (7.4-15.3) 

| RBTAW 7.5 5.9" 5.8 
: (7.3-7.7) (5.1-7.0) (4.4-6.2) 

: RMTAW 4.6 3.7 3.7 
: (3.8-4.7) (3.0-4.0) (2.9-4.3) 

| RBTAH 13.1 9.8" 12.1 
(12.7-14.1) (8.6-11.0) (11.5-12.9) 

i RMXTAH 14.4 10.0°* 12.7 
(13.5-14.6) (9.8-11.2) (11.0-14.2) 

RMTAH 14.2 8.67 11.8 
(13.5-14.3) (8.0-10.3) (11.0-13.3) 

| RELL 40.3 38.6" 43.7 
(39.8-41.2) (37.6-39.1) (42.6-44.6) 

| RHLL 442 414^ 44.8 
(43.9-45.4) (40.2-42.5) (42.3-48.1) 

P R2FL 6.1 6.0 5.5 
(5.6-6.8) (5.6-6.8) (4.5-6.8) 

R3FL 7.6 73 6.8 
(7.1-7.9) (5.9-8.4) (5.5-7.4) 

R3TL 9.6 8.6 8.3 
(9.6-11.1) (7.7-9.5) (7.0-9.2) 

R5TL 4.8 42 4.1 
(4.3-4.9) (3.7-4.7) (2.8-6.1) 
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average p-distance of 2% (range: 0—39^). 
From these genetic results, the Doi Chang 
samples are regarded as conspecific with 7. 
verrucosus or T. shanjing and heterospecific 
to the other Thai species. 
Morphological Study 

A. total of 20 adult males were used for 
the morphological comparisons (Table 2). 
For the SVL, the Doi Chang samples (T. 
verrucosus) were smaller than 7. uyenoi, but 
nearly similar to T. anguliceps. For the BW, 
the Doi Chang samples were the lightest. In 
character ratios, the Doi Chang samples had 
the smallest values for RHW, RBTAH, 





1.0 


PC2 (24.9295 of variance) 


-1.5 -1.0 5 .Ü 


PCI (27.67% of variance) 


MXTAH, RMTAH, RFLL, and RHLL, but, 
in contrast, the RAGD of the Doi Chang 
samples was the largest. For the RTRL, the 
Doi Chang samples were smaller than T. 
uyenoi but larger than T. anguliceps. 

The overall morphological differences in 
males were examined using PCA for three 
T. anguliceps, seven Doi Chang samples 
and 10 T. uyenoi. The first two principal 
components accounted for 52.59% of the 
total variation. The two-dimensional plot of 
PC2 versus PCI showed that all three 
species were completely separated from one 
another (Fig. 3). 


e 7. anguliceps 
B 7. uyenoi 
A J. verrucosus 


5 1.0 1.5 


FIGURE 3. Principal component analysis of Tylototriton showing the two dimensional plots of the first 
principal component (PC1) versus the second component (PC2). 
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TABLE 3. Morphological differences in adult males of the three Tylototriton species from this study. 


Character T. verrucosus 


Snout Truncate 


Cranial crest Narrow, weakly curved 


medially at posterior 


T. uyenoi 
Rounded to blunt 


Wide, weakly or rather 
curved medially at 


T. anguliceps 
Truncate 


Narrow, distinctly curved 
medially at posterior end 


end posterior end 
Surface of cranial Rather smooth Rough Rough 
crest 
Sagittal ridge on head Rather distinct Rather distinct Prominent 
Vertebral ridge Narrow, segmented Wide, segmented Narrow, weakly segmented 
Rib nodules Isolated but few Isolated, rounded, 12-16 Isolated, rounded, 14—15 


individuals connected, 
rounded but few 
individuals oval, 12-16 


Glandular skin Small granules, sparse 


Adpressed limbs Overlapping 
Brown to dark brown 


Color of ventral trunk 


Dull orange-brown to 
yellowish-brown 


Color of digit tips 


The three species had a rather similar 
coloration, with a brown to dark-brown 
background, although the Doi Chang 
samples tended to be the most blackish. The 
dorsal and ventral head, parotoids, vertebral 
ridge, rib nodules, limbs, vent, and tail were 
dull orange-brown to  yellowish-brown, 
although T. uyenoi had darker orange 
markings than the other species. The ventral 
tail ridge was brighter in coloration among 
these three species. However, some parts of 
the bodies had different colorations (see 
Table 3). The morphological differences 
between the three species are summarized in 
Table 3. 

When the Doi Chang samples (Fig. 4) 
were compared with 7. verrucosus, using 








Large granules, dense 


Overlapping 


Brown to dark brown 


Dark brown to black 


Large granules, dense 
Overlapping 


Partly dull orange-brown to 
yellowish-brown 


Dull orange-brown to 
yellowish-brown 


data taken from the  neotype from 
Longchuan County, Yunnan Province, 
China, as described by Nussbaum et al. 
(1995), they both showed similar color 
patterns of a mostly blackish background 
but orange to reddish-brown markings on 
the head, trunk, limbs, vent, and ventral 
rdge of the tail. They had a distinctly 
segmented vertebral ridge. The rib nodules 
were round and distinct However, the 
number of rib nodules in the Doi Chang 
samples ranged from 12-16, whereas in 7. 
verrucosus (Nussbaum et al, 1995) the 
nodules amounted to 13. Both taxa had 
similar small and sparse granules that were 
denser on the dorsal side than the ventral 
side. 
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These morphological results conformed 
to the genetic ones mentioned above, where 
the Doi Chang samples are regarded as 
conspecific with T. verrucosus or T. 
shanjing and heterospecific to the other Thai 
species. Therefore, we provisionally 
identified the Doi Chang samples as 
Tylototriton cf. verrucosus at the moment. 
Habitat Use 

To date, in Thailand, Tylototriton cf. 
verrucosus is known only from a single 
locality at Doi Chang, Chiang Rai Province, 
which is located on the Khun Tan mountain 
range (Fig. 1). All adult males of 
Tylototriton cf. verrucosus were found at 
night between 21.00—22.00 h in the pond 
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FIGURE 4. Male Tylototriton verrucosus from Doi Chang, Chiang Rai Province, Thailand. 





during the breeding season (Fig. 5). The 
bottom of the pond was covered with small 
pieces of debris; such as leaf litter. 
Surrounding areas were hill evergreen 
forest. The pond size was approximately 
400 cm long, 290 cm wide, and 33 cm in 
maximum depth. At 21.34 h on the 2™ June 
2019, the water and air temperatures were 
20.4 °C and 22.4 °C, respectively, with a 
relative air humidity of 88%. The water 
quality parameters were: pH 6.75; dissolved 
oxygen 1.48 mg/L; conductivity 375 uS/cm; 
and total dissolved solid 164 mg/L. The 
altitude was 1500 m amsl. No fish were 
observed. 
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FIGURE 5. Pond habitat of Tylototriton verrucosus at Doi Chang, Chiang Rai Province, Thailand. 


DISCUSSION 


The p-distance between species of 
Tylototriton was previously estimated at > 
3.3% (Nishikawa et al., 2013a), which was 
confirmed in the present study except for 7. 
shanorum and T. ngarsuensis with a smaller 
distance of 2%, the same as the intraspecific 
variation in T. uyenoi. Thus, the variation 
within the three clades (Clades A, B, and C) 
in this study could be judged as being at the 
intraspecific level and the taxonomic 
relationship between T. shanorum and T. 
ngarsuensis should be reexamined in a 
future study. 

Clade A included topotypic samples of 
T. verrucosus and T. shanjing, and we 
tentatively treat this clade as 7. verrucosus. 
Anderson (1871) did not specify the precise 
type locality, but only mentioned some areas 
located at Nantin, Momien, and Hotha 
valleys, in western Yunnan Province 
(China). Moreover, Anderson (1871) neither 


classified the type series nor mentioned 
where the described specimens were 
deposited. Nussbaum et al. (1995) 
designated a neotype of the species from 
Gongwa Village, Longchuan County, 
Yunnan Province, China, because the 
presumed Anderson’s syntype was 
considered as T. shanjing and did not 
correspond to A Anderson's original 
description of the dark brown-patterned 
specimen. Tylototriton verrucosus and T. 
shanjing occur in parapatry in Yunnan, but 
they share mitochondrial genomes (Zhang et 
al., 2007; Phimmachak et al., 20152). Thus, 
further intensive study is needed to clarify 
their actual species identities. However, 
from the morphological and molecular data, 
we tentatively treat 7. verrucosus and T. 
shanjing as the same species and use the 
former species name here, and so identify 
the Tylototriton species from Doi Chang, 
Chiang Rai Province, Thailand as T. 
Verrucosus. 
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Clade B could be treated as 7. anguliceps 
and Clade C as T. uyenoi without any doubt 
from the present results. Accordingly, we 
now confirm that four species of 
Tylototriton exist in Thailand. Tylototriton 
anguliceps, T. uyenoi and T. verrucosus are 
found in the Northern part, while 7. panhai 
is only found in the North-eastern part, and 
they are all geographically separated. 
Although 7. verrucosus 1s distributed on the 
Khun Tan mountain range, Chiang Rai 
Province, close to the range of T. 
anguliceps, until now no sympatry between 
the two species has been observed. Since the 
Tylototriton species use shallow and 
isolated ponds on mountain ranges at 
elevations > 1000 m ams! for breeding, each 
population may be isolated owing to the 
presence of lowlands between the mountain 
ranges where they occur. 

During the breeding season, adult males 
of T. verrucosus from Doi Chang use 
similar microhabitats in ponds at more than 
1000 m amsl, as per previous reports in 7. 
anguliceps (Le et al, 2015), T. uyenoi 
(Pomchote et al., 2008: as the Type I form 
of T. verrucosus, Dowwiangkan et al., 
2018), and 7. panhai (Pomchote et al., 
2008; Khonsue et al., 2010: as the Type II 
form of T. verrucosus). 

We did not find any female newts in the 
three sampled locations. Female newts may 
enter breeding ponds later or they may stay 
in the water for only a short period and then 
return to land (Grismer et al., 2018). Males 
were found in a greater number than females 
in T. uyenoi (Dowwiangkan et al., 2018) and 
in some populations of T. himalayanus 
(Wangyal and Gurung, 2012: as T7. 
verrucosus). However, an equal sex ratio 
was observed in T. podichthys (Phimmachak 
et al., 2015b). Thus, a longitudinal study is 
needed for a population study in the Thai 
Tylototriton species. 


In the amphibian wildlife protection list 
of Thailand, the only protected newt species 
is T. verrucosus (Department of National 
Parks, Wildlife and Plant Conservation). 
The other three candidates (7. uyenoi, T. 
anguliceps, and T. panhai) are not protected. 
Moreover, the 7. verrucosus population, 
which is only found at Doi Chang, Chiang 
Rai Province, may be in risk of decline 
because the Doi Chang area is not protected 
by law. Thus, we suggest that this species 1s 
threatened because of habitat alteration and 
human activities, such as modifying the area 
into tourist attractions, tourist 
accommodation, and commercial crop 
cultivation plots, like Doi Chaang Coffee. 

Our study confirms that T. verrucosus 
sensu stricto 1s still 1n existence, even after 
recent and drastic taxonomic changes, in 
Thailand and is expected to be found in 
other areas of Thailand. However, their 
habitats need to be protected by future 
conservation plans. 

The experimental protocol was approved 
by the Animal Care and Use Committee of 
Faculty of Science, Chulalongkorn 
University (Protocol Review No. 1723008 
Renew). This research was permitted by the 
National Research Council of Thailand 
(NRCT) and the Department of National 
Parks, Wildlife and Plant Conservation 
(DNP). 
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